Power Spectrum of the density of cold atomic gas 
in the Galaxy towards Cas A and Cygnus A 
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ABSTRACT 

We have obtained the power spectral description of the density and 
opacity fluctuations of the cold HI gas in the Galaxy towards Cas A, 
and Cygnus A. We have employed a method of deconvolution, based on 
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CLEAN, to estimate the true power spectrum of optical depth of cold HI 
gas from the observed distribution, taking into account the finite extent 
of the background source and the incomplete sampling of optical depth 
over the extent of the source. We investigate the nature of the underlying 
spectrum of density fluctuations in the cold HI gas which would be 
consistent with that of the observed HI optical depth fluctuations. These 
power spectra for the Perseus arm towards Cas A, and for the Outer arm 
towards Cygnus A have a slope of 2.75±0.25 (3a error). The slope in 
the case of the Local arm towards Cygnus A is 2.5, and is significantly 
shallower in comparison. The linear scales probed here range from 0.01 
to 3 pc. We discuss the implications of our results, the non-Kolmogorov 
nature of the spectrum, and the observed HI opacity variations on small 
transverse scales. 

Subject headings: radio lines: ISM - methods: data analysis - tech- 
niques: image processing, interferometric - interstellar medium: struc- 
ture - turbulence 
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1. INTRODUCTION 

The cold atomic gas in the Galaxy has been 
observed for over forty years in the 21 cm HI 
line. These early single dish and interferomet- 
ric observations established the presence of 
cold atomic gas in the Galaxy possibly in the 
form of discrete 'clouds' in pressure equilib- 
rium with the warm neutral medium (Clark, 
Radhakrishnan, & Wilson 1962, Clark 1965, 
Radhakrishnan et al. 1972, Field, Goldsmith, 
& Habing 1969, Field 1973). These HI clouds 
with an estimated mean spin temperature of 
~ 80 K, and a mean density of ~ 30 cm -3 
were thought to be a few parsec in size. This 
simple picture of the cold gas has changed 
significantly during the last two decades due 
to high resolution 21 cm line absorption mea- 
surements of the Galaxy made using aperture 
synthesis telescopes (Greisen 1973, Lockhart, 
& Goss 1978, Bregman et al. 1983, Schwarz 
et al. 1986, 1995). Fine structures in the 
atomic hydrogen gas have been observed over 
a wide range in length scales down to about 
0.1 pc (Crovisier et al. 1985, Kalberla et al. 
1985, Bieging et al. 1991, Roberts et al. 1993, 
Schwarz et al. 1995). Some lines of sight 
which have also been explored through VLBI 
and multi-epoch pulsar observations have re- 
vealed significant opacity variations in the 
atomic gas on transverse scales down to ~ 
10 AU (Dieter et al. 1976, Diamond et al. 
1989, Deshpande et al. 1992, Frail et al. 1994, 
Davis et al. 1996, Heiles 1997, Faison et al. 
1998). The observed opacity variations had 
raised concerns regarding the origin, lifetime, 
and pressure equilibrium with the rest of the 
interstellar medium of small scale structure 
presumed to cause them. Related issues in 
this context are coherent length scales of the 
underlying structures, the role of turbulence 
in producing these, and the power contained 



in such structures as opposed to that in large 
scale structures. 

Some of these issues are best addressed 
through a power spectrum analysis of the HI 
density fluctuations at the relevant range of 
scales. For the warm HI gas observed in emis- 
sion, this is rather straightforward since the 
observed 21 cm line intensity is directly pro- 
portional to the column density of the HI 
gas. Thus, the observed line visibilities as a 
function of spatial frequency define the power 
spectrum of the HI density fluctuations. For 
scales of 50 to 200 pc this has been carried 
out by Green (1993), and a useful analytical 
formulation related to such measurements is 
discussed by Lazarian (1995). The measure- 
ment of HI absorption distribution, on the 
other hand, can not be directly interpreted in 
terms of the power spectrum of the HI density 
fluctuations, as the HI absorption depends on 
the integral along the sight line of the ratio of 
the number density to the spin temperature 
of the HI gas. Thus, the translation of the 
observed optical depth fluctuations to those 
in HI density necessarily involves an assump- 
tion about the dependence of the density on 
temperature (e.g., pressure equilibrium would 
imply constancy of the product of density and 
temperature ). Direct interpretations of the 
HI absorption measurements are made diffi- 
cult by the effects due to the finite extent and 
the shape of the background source. 

More recently, the opacity variations in 
the cold atomic gas have been imaged us- 
ing bright supernova remnants as background 
continuum sources (Beiging, Goss, & Wilcots 
1991, Roberts et al 1993, Schwarz et al 1995). 
Some of these observations provide a unique 
opportunity to study the opacity distribu- 
tions over a large range of spatial scales. In 
this paper, we attempt to derive an under- 
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lying power spectrum from the HI opacity 
distribution observed in the direction of two 
bright radio sources, Cas A, and Cygnus A. 
The power spectrum of the opacity distri- 
bution over about two orders of magnitude 
in spatial scales (as from the Cas A data) 
is found to follow a power law. The opti- 
cal depth images used in the current analysis 
are discussed in Section 2. The power spec- 
trum analysis and application of CLEAN for 
high-dynamic-range spectral estimation are 
discussed in Section 3. This CLEAN proce- 
dure, to recover the true spectrum from the 
observed spectrum, takes into account the ex- 
tent and shape of the background source. In 
section 4, the estimated power spectra of the 
cold HI gas in the Galaxy towards Cas A are 
presented. In this section, we also present 
the structure function of the observed opti- 
cal depth distribution in the Galaxy towards 
the direction of both the sources. The pos- 
sible relation between the opacity spectrum 
and the underlying spectrum of HI density 
fluctuations is examined in Section 4.1. The 
implications of the derived power spectra of 
the cold gas for the interstellar medium are 
discussed in Section 5. 

2. OPTICAL DEPTH IMAGES 

The distribution of cold HI gas in the 
Galaxy towards the supernova remnant Cas- 
siopeia A (In ~ 111°, hn ~ —2°) was ob- 
served by Bieging, Goss, & Wilcots (1991; 
hereafter BGW) using the Very Large Array 
(VLA) of the National Radio Astronomy Ob- 
servatory. They present optical depth images 
at an angular resolution of 7" and a spectral 
resolution of 0.6 km s -1 . Cas A is at a dis- 
tance of ~ 3 kpc from the Sun and is be- 
lieved to be in, or, beyond the Perseus spiral 
arm. The optical-depth image cube presented 



by BGW covers a velocity range from ~ -30 
km s _1 to -56 km s _1 (Vj ar ). The HI ab- 
sorption in this velocity range is caused by 
the interstellar gas in the Perseus spiral arm. 
An optical depth image in one of the veloc- 
ity channels (at V; sr = -38.8 km s _1 ) from 
this cube is reproduced in Fig. 1. The mor- 
phology of the absorbing gas at an adopted 
distance of ~ 2 kpc to the Perseus arm (spa- 
tial resolution ~ 0.07 pc), is dominated by 
filaments, loops, and arcs. Recently, we used 
the VLA to observe the distribution of cold 
gas in the Galaxy toward the direction of the 
powerful extragalactic radio source Cygnus A 
(In ~ 75°, bn ~ 5°). The observations were 
carried out with the VLA in the B configu- 
ration (maximum baseline ~ 11 km) on De- 
cember 26, and 27, 1999. A bandwidth of 
1.56 MHz was used with 256 channels, and 
Hanning smoothing, resulting in a spectral 
resolution of 1.29 km s _1 . The synthesized 
beam was 3". 8 x 3". 6 at a position angle 
of -2°. 8. Bandpass calibration was achieved 
by frequency switching on Cygnus A. The 
channels from the flat portion of the spec- 
trum and free from line absorption were used 
to estimate continuum visibilities. Spectral- 
line visibilities in each of the relevant spec- 
tral channels were obtained by removing the 
estimated continuum contribution to the vis- 
ibilities. The images constituting a spectral 
image cube, as well as the continuum image 
were constructed and CLEANed. The rms 
noise in each of the spectral channels and in 
the continuum image was 0.1 Jy/beam, lim- 
ited by dynamic range. The optical depths 
(r) were estimated in the standard way at all 
positions where the line intensity and the con- 
tinuum intensity were greater than five times 
the rms value in the respective images. Some 
examples of optical depth images produced 
this way are shown in Figs. 2, and 3. The 
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Fig. 1. — The HI optical depth image at a V; sr = —38.8 km s _1 towards Cas A from BGW. This 
VLA image has an angular resolution of 7" with a sampling interval of 1".6. The optical depth 
grey scale ranges from 0.5(white) to 5 (black). The contours are in units of 0.1 from 0.1 to 1, and 
in units of 1 from 1 to 5 in optical depth. The synthesized beam is shown in the bottom right 
hand corner. The velocity resolution of this image is 0.6 km s _1 . 
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Fig. 2. — The HI optical depth image of the Local arm towards Cygnus A produced from the VLA 
observations described in this paper. The top right hand corner indicates the Vi sr - The bottom 
left hand corner of the first panel shows the synthesized beam. The optical depth grey scale ranges 
from (white) to 0.4 (black). The contours are in units of 0.01 from 0.01 to 0.1, and in units of 
0.05 from 0.1 to 0.5 in optical depth. The angular resolution of ~ 3. "7 corresponds to a linear 
resolution ~ 0.01 pc for an assumed distance of 0.5 kpc to the Local arm. The sampling interval 
in these images is 1", and the velocity resolution is 1.29 km s _1 . 
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Fig. 3. — The HI optical depth image of the Outer arm towards Cygnus A obtained from the 
VLA. The contours are as in Fig. 2. The angular resolution of ~ 3. "7 corresponds to a linear size 
~ 0.27 pc for an assumed distance of 15 kpc to the Outer arm. 
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Fig. 4. — The upper and lower frames show the HI optical depth spectra towards the northwest and 
the southeast peaks of Cygnus A respectively. Note the large optical depth values corresponding 
to the Local (V; sr ~ km s -1 ), and the Outer (Vi sr ~ -85 km s" 1 ) arms. In contrast, the Perseus 
arm shows a weak HI absorption (r < 0.05) at Vi sr ~ -50 km s _1 . 
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double-lobed morphology is simply a reflec- 
tion of the radio continuum morphology of 
Cygnus A. The optical depth images in Figs. 

2, and 3 provide examples of the distribution 
of cold atomic gas in the Local, and the Outer 
spiral arms, respectively. The spatial resolu- 
tions of these images correspond to ~ 0.01 
pc, and 0.3 pc at the Local, and the Outer 
arms respectively. The assumed distances to 
the Local, and the Outer arms are 0.5, and 
15 kpc, respectively. There is very little HI 
absorption (peak optical depth < 0.05) in 
the Perseus arm towards Cygnus A (Fig. 4). 
Consequently, except towards the two bright- 
est positions in the two lobes of Cygnus A, 
the optical depth of the Perseus arm remains 
unmeasured towards most of the area cov- 
ered by the two lobes. These images are not 
displayed. For the same reason, we also do 
not use the Perseus arm optical depth images 
towards Cygnus A in the structure function 
analysis discussed in Section 4. 

3. POWER SPECTRUM ANALYSIS 

For the first time, the optical depth im- 
ages in the Cas A direction provide possibil- 
ities of studying the distribution of cold HI 
over a usefully wide range of spatial scales si- 
multaneously. Given the extent and the reso- 
lution of these images, a detailed estimate of 
the power spectrum of the opacity distribu- 
tion over nearly two decades in spatial scales 
is possible. In this section, we attempt such 
an estimate and present details of our analy- 
sis. 

In an ideal situation where a fully sam- 
pled opacity distribution is available over an 
extent of 9 max and with a sampling interval 
A9, it is possible to simply Fourier trans- 
form the distribution and estimate the power 
spectrum over a frequency range l/9 max to 



1/(2 A0). In reality, the opacity distribution 
is not fully sampled (see Fig. 1), making it 
difficult to estimate the true power-spectrum 
directly. As is clear from the distribution in 
Fig. 1, the incomplete sampling is primar- 
ily due to the finite extent of the background 
source, Cas A. In principle, a smaller sec- 
tion of the image where sampling is complete 
can be chosen; however, reduced resolution 
in the spectrum then results. Unfortunately, 
some incompleteness exists in most of the op- 
tical depth images within the source extent 
where the continuum background or the line 
intensity is too weak to allow a reliable es- 
timation of the optical depth. If the optical 
depth at the unsampled locations in the im- 
age is assigned 'zero' value, the Fourier trans- 
form of the image gives a spectrum that is 
a convolution of the true spectrum with the 
Fourier transform of the measurement mask 
or the sampling spectrum. The measurement 
mask is defined by ones and zeros at the lo- 
cations that are sampled, and unsampled, re- 
spectively. This convolution can significantly 
modify the spectrum, thus making the nature 
of the true spectrum less apparent. Some 
form of deconvolution operation is therefore 
required to recover the true spectrum, partic- 
ularly when high dynamic range estimation is 
a primary consideration (as would be while 
studying steep power-law spectra). 

3.1. Application of CLEAN for high- 
dynamic-range spectral estimation 

The situation in obtaining the 'true' spec- 
trum from the modified or dirty spectrum is 
analogous to that of obtaining a clean image 
from a dirty image routinely encountered in 
synthesis imaging (Ryle 1960). There, the 
'dirty image' obtained from a telescope is 
the convolution of the 'clean image' with the 
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'point spread function' (PSF) of the telescope 
and the dirty image needs to be deconvolved 
with the PSF in order to obtain the true (or 
clean) image. The deconvolution attains even 
more importance when the telescope is an un- 
filled aperture. In such a situation a technique 
called CLEAN (Hogbom 1974) is commonly 
used to obtain high-dynamic-range images. 
In the present situation, the dirty spectrum is 
to be deconvolved using the spectrum of the 
mask in order to obtain the true spectrum. 
The application of CLEAN in the power spec- 
tral analysis of non-uniformly sampled pul- 
sar timing data was developed by Deshpande 
et al. (1996). In that case a 1-dimensional 
CLEAN was developed. In the present con- 
text a 2-dimensional variant of this applica- 
tion is used. Details of the application of 
CLEAN in the context of power spectral anal- 
ysis can be obtained from Deshpande et al. 
(1996). Only a brief summary of the proce- 
dure used in the present context will be given 
here. 

The 2-dimensional r image (for e.g., Fig. 
1) which is 223 x 223 pixels is padded with 
zeros to obtain a 512 x 512 image. A mask 
of the same dimensions is also prepared. This 
padding helps in improving the sampling by 
a factor of ~2 in the spectral domain which 
improves the performance of CLEAN quali- 
tatively, although at the expensive additional 
computation. Both the image and the mask 
in the image domain are Fourier transformed 
to obtain the image spectrum and the mask 
spectrum in the spatial-frequency domain, re- 
spectively. The spectra so obtained are her- 
mitian symmetric (as are visibilities and un- 
like the images in the case of synthesis imag- 
ing) and the CLEAN algorithm was modi- 
fied to account for this modification: (i) Dur- 
ing each iteration of CLEAN, while search- 



ing for the maximum in the spectrum, only 
the amplitudes were considered. But, while 
removing CLEAN components, the complex 
nature (amplitude, and phase) is taken into 
account, and, (ii) the search for the maxi- 
mum was made only over one half of the im- 
age spectrum. Subsequently, subtraction of a 
scaled version of the mask spectrum was per- 
formed over both halves of the image spec- 
trum after accounting for the hermitian sym- 
metric contribution. This procedure ensures 
the hermitian symmetry in the CLEAN spec- 
trum. The CLEAN complex components are 
restored with a half-cycle cosine of half power 
width close to the original spectral resolution 
which is added to the residuals. The square 
of the amplitude of the resultant CLEAN 
complex spectrum is then used to estimate 
the power spectrum. Finally, an azimuthally- 
averaged 1-d power spectrum is obtained from 
such a 2-d CLEAN power spectrum to im- 
prove the spectral definition; of course any 
information on anisotropy is lost. This proce- 
dure was applied to a series of optical-depth 
images corresponding to individual velocity 
channels as well as those integrated over wider 
velocity ranges. Given the continuous nature 
of the spectra, one may legitimately ques- 
tion the applicability of the CLEAN decon- 
volution procedure (Schwarz 1978), since the 
CLEAN assumes a generally 'empty' domain 
with only well separated narrow spectral com- 
ponents in the 'true' spectrum. Interestingly, 
the often encountered steep power-law spec- 
tra (power oc q~ a , where q is the spatial fre- 
quency) are 'red' in nature, i.e., a is posi- 
tive. When viewed on a linear-linear scale, 
such spectra would be highly peaked around 
'zero' frequency and relatively 'empty' else- 
where, providing a favorable situation for the 
CLEAN. 
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The azimuthally-averaged mask spectrum 
corresponding to the image in Fig. 1 is a 
featureless power-law (power oc q _a , where 
q is the spatial frequency) with a slope (or 
power-law index) of 2.9 (—a). This spec- 
trum can be understood by noting that the 
optical-depth image in Fig. 1 defines an ap- 
proximately circular area (the morphology of 
the background source Cas A) as the sampled 
distribution. The Fourier transform of a uni- 
formly weighted perfect circular patch is an 
Airy pattern, i.e., a Bessel function of I order 
divided by q. For large values of its argument 
(q), the square of the Airy pattern decreases 
as q~ 3 . Hence, the observed power-law slope 
of the mask-spectrum is not surprising. As an 
example of the improvement brought about 
by the CLEAN procedure, Fig. 5 shows the 
azimuthally-averaged versions of the dirty- 
spectrum (i.e. pre-CLEAN, shown as dashed 
line) and the clean-spectrum (i.e. CLEANed 
and restored, shown as solid line) for the im- 
ages similar to that in Fig. 1. In both cases 
the power spectra from eleven channels in the 
velocity range -41.4 km s" 1 to -34.9 km s _1 
were averaged to improve the spectral defi- 
nition, since the individual noisy spectra are 
similar. There are two distinct differences be- 
tween the two spectra. First, the true spec- 
trum has a knee at log(q/q ma:r ) ~ —0.5 corre- 
sponding to the size of the synthesized beam 
~ 7". The power spectrum drops sharply be- 
yond these frequencies. This 'bowl' is filled in 
by the 'leakage' of power from lower frequen- 
cies in the dirty spectrum. Second, the slope 
of the true spectrum is 2.75 while that of the 
dirty spectrum is 2.4 (both the slopes calcu- 
lated between -0.5 and -1.5 in log(q/q max )). 
The dirty spectrum is expected to be flatter 
than the CLEAN spectrum since the latter 
has been convolved with the mask spectrum 
leading to a mixing (spread) of power from 



smaller to larger spatial-frequencies. 

It may be noted that although the process- 
ing described so far involves 2-dimensional 
images/spectra, the results are presented as 
azimuthally averaged versions. This averag- 
ing improves the definition of the spectra in 
a monotonically increasing fashion towards 
higher spatial frequencies consistent with the 
increasing number of independent samples 
available for the averaging. Thus the uncer- 
tainties are relatively high towards the lower 
spatial frequency end of the spectrum, which 
is therefore excluded from any quantitative 
estimation (e.g. slope). At the other end, the 
shape of the 'bowl' is a result of the steepening 
of the spectrum due to beam smoothing and 
a flattening of the spectrum as the contribu- 
tion due to the uncertainties in the estimation 
of optical depth becomes significant. Consid- 
eration of these two effects defines the range 
of spatial frequencies over which the spectral 
characteristics can be estimated reliably. The 
CLEAN spectrum viewed over this range is 
describable by a simple power law of the na- 
ture q~ a , where q is the spatial frequency, and 
a is the power-law index. 

4. POWER SPECTRA FROM CAS 
A AND CYGNUS A 

The spatial power spectra of optical-depth 
distribution towards Cas A are shown in Fig. 
6. These spectra represent average descrip- 
tion over the velocity intervals -51.6 to - 
45.2 km s _1 and -41.4 to -34.9 km s^ 1 sep- 
arately. These velocity ranges correspond to 
the Perseus arm and arise from a distance of 
2 kpc. In Fig. 6 there are three spectra and 
are obtained as follows : (1) In the first case 
(Case A), a power spectrum is obtained for 
each of the 11 velocity channels in the cho- 
sen velocity range (for e.g., -51.6 to -45.2 km 
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Fig. 5. — Demonstration of the CLEAN procedure in power spectra. The x-ordinate is logarithm 
of the spatial frequency in units of the limiting frequency (q mox ). The limiting frequency is half of 
the sampling frequency. The sampling interval in the Cas A optical depth images used here is l."6. 
The top curve (dashed line) is a dirty power spectrum obtained after averaging the power spectra 
of 11 channels in the velocity range Vi sr = -41.4 km s _1 to -34.9 km s -1 . The bottom curve (solid 
line) is the true power spectrum obtained after CLEANing. Note the appearance of the 'knee' in 
this spectrum at log(q/q maa; ) ~-0.5. This effect is due to the finite size (7") of the synthesized 
beam. The slopes of the dirty and true spectra are 2.4 and 2.75 respectively, estimated over the 
x-range -0.5 to -1.5. 
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s -1 ). The resultant 11 power spectra were 
found to be similar in slope and differing in 
the spectral density consistent with the opti- 
cal depth variations in the images used. They 
were combined with uniform weighting giv- 
ing an average description of opacity distri- 
bution viewed over a narrow channel (width 
0.6 km s _1 in the present case). In Fig. 6 
there are two spectra obtained in this fash- 
ion corresponding to the two velocity ranges 
-51.6 to -45.2 km s _1 (indicated by dots and 
dashes) and -41.4 to -34.9 km s -1 (solid line). 
These two spectra are less noisy (owing to the 
averaging over 11 channels) and follow each 
other well in the log(q/q max ) range of -1.5 
to -0.5. It is not surprising that the agree- 
ment of the power spectra in these two veloc- 
ity ranges is excellent given the same range 
of spatial scales covered. In the second case 
(Case B), the optical depth images over 11 
velocity channels of a chosen velocity range 
(-41.4 to -34.9 km s _1 ) were first averaged 
to obtain one image for which a power spec- 
trum was estimated. A given pixel in the av- 
erage image was blanked when the average 
could not be performed over all the 11 veloc- 
ity channels (i.e., when one or more of the 
channel images had the corresponding pixel 
blanked). This reduces, however, the total 
number of 'valid' pixels in the resultant im- 
age. This effect was minimized for the veloc- 
ity range -41.4 to -34.9 km s -1 . This spec- 
trum is shown in Fig. 6 CIS db dashed-line and 
is easily recognized as the noisiest of the three 
spectra shown. The slopes estimated over the 
l°g(q/ Qmax) range of -0.5 and -1.5 for all the 
three spectra are equal within the noise, with 
a power law slope of 2.75±0.25 (3a error). 
The average power in the Case B spectrum 
(dashed line) over this range in log(q/q maa; ) 
is about a factor of three higher than the two 
Case A spectra. The implications of this en- 



hanced power will be discussed in Section 5. 

A similar analysis was also attempted for 
the data in the direction of Cygnus A . How- 
ever, given the very limited sampling in this 
case, a structure function should provide a 
more useful description of the opacity distri- 
bution than an azimuthally averaged spec- 
trum even after CLEAN. Hence, we have ob- 
tained the structure function for the observed 
optical depth distribution towards Cygnus A 
and compared this with the Cas A results. 
The structure function of the optical depth 
is defined as D T (r) =< [t(x) — t(x — r)] 2 >, 
where t(x) is the optical depth at point x, and 
the estimation is based on only those points in 
the image where the optical depth estimates 
are considered to be reliable. It can be shown 
that for a distribution having a power law 
spectrum (oc q~ a ), the structure function will 
also be a power law (oc r 2/3 , where r oc l/q), 
such that a = 2+2f3 (Lee & Jokipii 1975). 

In Fig. 7, the square root of the struc- 
ture function (i.e., the rms fluctuation) of the 
optical depth distribution in the Perseus arm 
towards Cas A (indicated by xs) is compared 
with those in the Local (indicated by cir- 
cles) and Outer arms (indicated by squares) 
towards Cygnus A. The spatial scales indi- 
cated are appropriate to the assumed dis- 
tances of the absorbing HI gas. The effect 
of the beam smoothing is apparent towards 
the smaller spatial scales (over about half an 
order of magnitude) in each of the structure 
functions. Towards the larger spatial scales, 
the estimation suffers from an inherent limi- 
tation in having a limited number of indepen- 
dent realizations of the large scale variations 
that are required to obtain a true average esti- 
mate. The gap seen in the structure function 
corresponding to the Outer arm of Cygnus A 
results from inadequate statistics consistent 
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Fig. 6. — The spatial power spectra of optical depth distribution towards Cas A. The spectrum 
indicated by the solid line is obtained as in Case A (Section 4) and corresponds to the velocity 
range -41.4 to -34.9 km s _1 . The spectrum indicated by the dot-dashed line also corresponds to 
Case A, but to the velocity range -51.6 to -45.2 km s _1 . The third spectrum (dashed line) is 
obtained as in Case B (see Section 4) and corresponds to the velocity range -41.4 to -34.9 km 
s _1 . All the three spectra have a slope of 2.75±0.25 (3a error) over the log(q/q max ) range of -0.5 
and -1.5. However, the Case B spectrum has three times higher power than the Case A spectrum 
over this range of log(q/q ma:r ). 
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with the absence of optical depth estimates 
over a large fraction of the north-west lobe 
(see Fig. 3). Even with these limitations 
the reliably estimated sections of the struc- 
ture functions indicate that the slopes of two 
of the three determinations (the Perseus arm, 
and the Outer arm) are consistent (f3 ~0.4). 
In addition, this value of /3 is consistent with 
the power law index (a = 2.75) for the gas in 
the Perseus arm (Fig. 6). The slope in the 
case of the Local arm is significantly shallower 
(P ~0.25). In addition, the fluctuations in the 
opacity (and its mean value) are an order of 
magnitude higher for the gas associated with 
the Perseus arm of Cas A as compared to the 
Outer arm HI in the direction of Cygnus A. 
The vertical heights from the Galactic plane 
of the HI gas observed in these two cases are 
different, viz., ~ 70 pc for the gas in the 
Perseus arm towards Cas A, and ~ 1300 pc 
for the gas in the Outer arm towards Cygnus 
A. However, it is unclear how much of the 
differences seen in the structure functions cor- 
responding to these arms are attributable to 
the differences in the heights from the Galac- 
tic plane of the absorbing HI gas in these two 
arms. 

4.1. Implications for the distribution 
of HI density and for pressure 
equilibrium 

The derived (CLEAN) power spectra shown 
in Fig. 6 represents p T (k), the spatial power 
spectrum of the optical depth distribution. 
The optical depth, r, is related to the col- 
umn density (N#) and (density weighted har- 
monic) mean spin temperature (T s ) of the HI 
gas as t oc N H /T S . The column density of 
HI is related to the volume density (n#) as 
Nh = n#l, where, 1 is the line of sight path 
length of the HI gas. The pressure of the HI 



gas (P) viewed locally is related to the volume 
density and kinetic temperature (T^) as Poc 
n^Tfc, assuming it to be a perfect gas. For the 
cold gas in the Galaxy, T s m T&. If the differ- 
ent density regions are in pressure equilibrium 
(i.e., nj/ oc 1/T&) with each other then this 
implies that r oc n 2 H . In general, if the equa- 
tion of state for the cold HI gas is given by 
Tfc ocn^, then the optical depth will be pro- 
portional to n^ 7 . Pressure equilibrium cor- 
responds to 7 = -1. It is not straightforward 
to predict the power law index of the den- 
sity spectrum based on that of the observed 
optical depth. However, it is possible to ex- 
plore the relation between these two through 
suitable Monte Carlo simulations involving an 
assumed power law spectrum of density and 
a value of 7. Our attempts to do so have in- 
dicated that the slope of the computed power 
spectrum of opacity is nearly the same as that 
of the assumed slope of the 3-d power spec- 
trum for the HI density. The two slopes differ 
marginally if in the assumed density distri- 
bution, the density fluctuations far exceeded 
the mean density. Otherwise, the resultant 
slope showed no dependence on the assumed 
value of 7, indicating that the assumption 
of pressure equilibrium is reasonable. This 
independence of the slope on the assumed 
value of 7 is not surprising since any non- 
linear dependence of r on nn implied by the 
equation of state (r oc n^ 7 ) reduces to the 
first-order to a dominant linear dependence 
(St oc (l-j)Sn H ) if the fluctuation in density 
(Sjih) is small compared to the mean density 
(<n H >). Hence, the nature of the observed 
spectrum of r distribution must closely de- 
scribe the corresponding spectrum of the HI 
density. 
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spatial scale (pc) 

Fig. 7. — The structure function towards Cas A and Cygnus A. The square root of the structure 
function (i.e., the r.m.s. value) of the optical depth distribution in the Perseus arm towards Cas A 
(indicated by xs) is compared with those towards Cygnus A in the Local (indicated by circles) and 
the Outer (indicated by squares) arms. The spatial scales correspond to the assumed distances of 
0.5, 2, and 15 kpc for the Local, the Perseus, and the Outer arms respectively. See Section 4 for 
further details. 



15 



5. DISCUSSION 

As can be seen in Fig. 6, the spectrum 
obtained in Case B (spectrum of the optical 
depth averaged over several velocity channels) 
has about a factor of 3 higher power than 
the Case A spectrum (average of the spectra 
of several velocity channels). This difference 
indicates that there is partial correlation in 
structures as a function of velocity. The spec- 
tra are normalized suitably such that if the 
structures were fully correlated across N chan- 
nels, we would observe an increase in power 
by a factor of N. Hence, the observed increase 
by a factor of three in power indicates a cor- 
relation scale of ~ 2 km s _1 in velocity. This 
correlation scale in velocity is expected since 
the typical value of turbulence for the cold HI 
gas in the Galaxy is ~ 4 km s _1 (Radhakrish- 
nan et al. 1972). 

The power spectrum of density fluctua- 
tions in the cold atomic gas towards Cas A 
has a slope of ~ 2.75. Adopting a distance of 
2 kpc for this absorbing gas, the linear scales 
probed range between 0.07 and 3 pc, corre- 
sponding to the resolution (7") and the size 
of Cas A (~5') respectively. Compared to 
these transverse scales, the longitudinal (or, 
the line of sight) scales corresponding to each 
of the channels (i.e., velocity resolution) must 
be significantly larger. Otherwise, the power 
spectrum obtained from the integrated im- 
age (Case B) should have a different (steeper) 
slope compared to that for a single channel 
image (Case A). However, as has been demon- 
strated in Fig. 6, there is no significant dif- 
ference in the two slopes. Thus, the opacity 
image in a single channel corresponds to a lon- 
gitudinally integrated version of the underly- 
ing 3-dimensional opacity distribution. Com- 
bining this result with the discussion in Sec- 
tion 4, the power spectrum represents the 3- 



d power spectrum of density fluctuations in 
the cold atomic gas over linear scales of 0.07 
to 3 pc. The HI emission analysis (Green 
1993) provided a similar slope, over the larger 
scales of 50 to 200 pc. Assuming no changes 
in the slope of the power spectrum between 
these two ranges, a constant slope is implied 
over 3 orders of magnitude in linear scales. 
It is noteworthy that the estimated slope of 
2.75 for the power spectrum differs signifi- 
cantly from the Kolmogorov slope (11/3 = 
3.67) expected for the 3-dimensional power 
spectrum of a turbulence-dominated incom- 
pressible fluid. 

In a recent paper, Lazarian, & Pogosyan 
(2000) draw attention to the modification of 
the HI density spectrum due to the veloc- 
ity field. They point out that the velocity 
range needs to be thick in order to obtain the 
true density power spectrum. They estimate 
a power law index close to the Kolmogorov 
value for the power spectrum of density using 
the HI emission data of the Galaxy (Green 
1993) and of the Small Megallanic Clouds 
(Stanimirovic et al. 1999). Being aware of 
this complication, we have obtained the power 
spectra before (Case A), and after (Case B) 
averaging in velocity (Section 4). Within the 
errors of estimation, the two slopes agree (Fig. 
6), indicating that the effect of averaging over 
velocity is not a substantial effect. 

It may not be surprising that the slope of 
the power spectrum of density fluctuations in 
the cold atomic gas in the Galaxy is different 
from the Kolmogorov slope considering that 
: (a) the HI gas in the interstellar medium 
is compressible, and, (b) the Kolmogorov in- 
dex is expected if the energy input into the 
ISM at some large scale cascades down to 
small scales through a turbulence mechanism. 
However, the exact nature and the scales at 
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which the energy is input, and its cascade to 
smaller scales are determined by a range of 
processes in the interstellar medium such as 
the supernovae, and stellar winds. Although 
a qualitative discussion of many such aspects 
is available (Lazarian 1995), the quantitative 
implications of these processes on the slope of 
the power spectrum is quite uncertain. 

Another important implication of the de- 
rived power spectral description is to the long- 
standing puzzle of the observed opacity vari- 
ations on small transverse scales down to tens 
of AU. These opacity variations have been in- 
terpreted in the past as evidence for wide- 
spread small-scale structure in the HI distri- 
bution with volume densities orders of magni- 
tude higher than those implied by the parsec- 
scale component. It has been recently empha- 
sized (Deshpande 2000) that there is a need 
to recognize the nature of the actual quan- 
tity measured at small scales; almost all scales 
contribute to the measured opacity differences 
at small scales. Furthermore, if the power 
spectrum with a slope of 2.75 estimated for 
the Cas A Perseus arm data is extrapolated 
to these small scales, optical depth variations 
with an r.m.s. of 0.1 can be expected on scales 
of ~ 100 AU (Deshpande 2000). The recent 
HI absorption measurements using the Very 
Long Baseline Array and the Very Large Ar- 
ray indicate that 5 out of the 7 sources ob- 
served show opacity variations consistent with 
this extrapolation (Faison et al. 1998, Fai- 
son & Goss, 2000). However, 3C138 & 3C147 
show optical depth fluctuations an order of 
magnitude larger over linear sizes of 10 - 20 
AU indicating real two dimensional structures 
in HI opacity but, note that the associated 
spectra appear much shallower than in the 

CclS .A. CclSG. 

Our investigations in two different direc- 



tions (Cas A, and Cygnus A) have shown both 
similarities and differences in the derived de- 
scriptions in terms of slopes, and powers as- 
sociated with the power spectra (and struc- 
ture functions) of the HI distribution. For 
a relatively complete picture to emerge it is 
desirable that such studies are carried out in 
different directions in the Galaxy. 

The analysis software used here was based 
largely on a package developed by Jayadev 
Rajagopal with one of us (AAD). We would 
like to thank Rajaram Nityananda, and V. 
Radhakrishnan for useful discussions, and com- 
ments on the manuscript. National Radio As- 
tronomy Observatory is a facility of the Na- 
tional Science Foundation operated under co- 
operative agreement by Associated Universi- 
ties, Inc. 
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